Entanglement Spectrum of a Disordered Topological Chern Insulator 
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We investigate the behavior of a topological Chern Insulator (CI) in the presence of disorder, 
with a focus on its entanglement spectrum (EtS) constructed from the ground state. For systems 
with symmetries, the EtS was shown to contain information about the topological universality 
class revealed by sorting the EtS against the conserved quantum numbers. In the absence of any 
symmetry, we demonstrate that statistical methods such as the level statistics of the EtS can be 
equally insightful, allowing us to distinguish when an insulator is in a topological or trivial phase and 
to map the boundary between the two phases. The phase diagram of a CI is computed as function 
of Fermi level (Ep) and disorder strength using the level statistics of the EtS and energy spectrum 
(EnS), together with a computation of the Chern number via an efficient real-space formula. 

PACS numbers: 63.22.-m, 87.10.-e,63.20.Pw 



Topological insulators (TI) are materials that do not 
conduct electricity in the bulk but display conducting 
edge channels. CIs represent a particular class of TIs [l] 
that have broken time-reversal symmetry. They haven't 
been observed yet experimentally, but a time-reversal in- 
variant version has been proposed [2]-[4] and discovered 
[5]. Since then, the field of TIs became increasingly pop- 
ular. The central claim of the field, and the basis for most 
potential applications, is the robustness of TIs' proper- 
ties against imperfections. Our Letter contributes to the 
ongoing research [6l? - HI] on disorder effects in TIs and 
gives two practicle tests to determine if the ground-state 
of a disordered CI is in the topological or trivial phase. 
The tests involve only the ground-state wavefunction. 

While our work concentrates on CIs, it addresses a 
broader question: Given the ground state of a Hamil- 
tonian, how much information can we extract about its 
"topological" universality class? Ref. [12] suggested that 
the answer resides in the entanglement spectrum^ i.e. the 
full set of eigenvalues of the reduced density matrix. For 
the i^=| Fractional Quantum Hall (FQH) states, the EtS 
levels and their multiplicities, when plotted versus the an- 
gular momentum, match the levels and multiplicities of 
the edge modes p!2HT4] . The EtS also captures the low- 
energy physics of gapless spin chains [15 and for topo- 
logical insulators it exhibits analogs of the physical edge 
state spectra ITBHTOl . 



mental symmetries of the Hamiltonians [20] . This line of 
thinking has had success in both many-body systems as 
well as in the theory of Anderson localization. We adopt 
it here and apply it to the EtS rather than to the EnS. For 
disordered CIs, we show that the EtS gives clear signa- 
tures of whether the CI is in the topologically nontrivial 
or in the simple Anderson insulator state. Such signa- 
tures, present in the ground state alone, are important 
for developing tools to attack the interacting many-body 
problem with disorder since statistical analysis of the EnS 
is impossible for large systems - diagonalization proce- 
dures only give a few low- lying energy states. The EtS 
of the ground-state, however, contains a large number of 
eigenvalues on which level statistics can be performed. 
We compare our results for the EtS with computations 
of the Chern number via a real-space formula and with 
results from conventional level statistics of EnS. 

We consider 2D lattice models with a = 1, . . . , i^ quan- 
tum states \x^a) per each site x and Hamiltonians: 
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Prior studies on EtS treated systems with translational 
invariance and the EtS of FQH states, spin chains, and 
topological insulators were plotted versus the momentum 
parallel to the cut. This is not generic and until now it 
was unclear if the EtS is useful when no symmetry is 
present. Without translational symmetry, what remains 
that is fundamental? Seminal thinking by Wigner gave 
us the answer: adopt a statistical view of spectra. Subse- 
quent work on random matrix theory revealed universal 
spectral properties that are dependent only on the funda- 



where the first term is a translationally invariant insu- 
lating Hamiltonian Hq and the second is a disorder po- 
tential V^. In 2D, Hq can display topological properties, 
manifested in the emergence of chiral edge modes. The 
number of stable chiral edge modes is equal to the C 
number of the occupied bulk states [3T]. By definition, 
a CI is an insulator with C^O. Its hulk states display 
a spectacular behavior, manifested in the persistence of 
extended states even when the disorder is on. The in- 
teresting physics of the CIs is due to these states - the 
edge modes are nothing but the delocalized bulk states 
terminating at the edge. 

In our calculations, we use the spin-up component of 
the Kane-Mele Hamiltonian [2j with A^^ = and Xso = 
T] — it (in order to connect with previous studies [10]): 
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FIG. 1. EnS and level statistics for a disordered (A) CI (C = 0.3i) and (B) trivial insulator (^ = 0.3), at disorder strengths 
W = (a) 3, (b) 5, (c) 8, (d) 11. Panel (A) also shows the phase diagram of the disordered CI inferred from panels (a)-(d), and a 
few histograms of the level spacings recorded at the indicated energies. The histograms are compared with Pgue(s) (blue line) 
and Ppoisson(s) (red lines) distributions. The blue lines overlaying the EnS are the variances of the energy-spacing distributions 
plotted on a scale between and 1. The dashed blue line represents the variance (^ 0.178) of Pque- 



where (x = ^o^x{t -\- iT]) with ax= the iso-spin of the site 
and cc, y are sites of a honeycomb lattice. This Hq dis- 
plays a topological phase for |?7| > \t\ tan |. For disorder 
we use uniform random entries uJx G [— ^, ^ • 

We first use the traditional level statistics analysis 
of the EnS to probe the extended/localized character 
of the bulk states. Fig. [l] shows the EnS of H^ at 
disorder strengths VK=3, 5, 8 and 11, when Hq is in 
topological (Fig lA) and trivial (Fig IB) phases. The 
energy levels are shown on the vertical axis for 10^ 
disorder configurations offset horizontally. Overlayed 
is the variance of the energy level spacings, collected 
at all energies using a small window. From the his- 
tograms, we see two regions where the level- spacing dis- 
tribution perfectly matches the Wigner-Dyson distribu- 
tion PG\jE{s) = ^^e~^^ '^ . The level-spacing variance 
(5^) — (s)^ at these energies converges to the variance of 
Pque— 0.178. We infer that these regions of level repul- 
sion contain extended states. In the rest of the spec- 
trum, the histograms match the Poisson distribution 
Pp(5)=e~* and the variance takes large values (0(1)). 
We infer that in these regions the states are localized. 
Upon increasing VK, the regions of delocalized spectrum 
(where the variance is exactly 0.178) converge towards 
each other to eventually collide and disappear. This is 
consistent with the "levitation and annihilation" phe- 
nomenon [To] and suggests the phase diagram shown in 
Fig. [T]A. In contrast, if we start from the normal insula- 
tor and increase W, all the states localize (see Fig. lip) 
at any disorder strength and there is no diffusive regime 
(a small variance in Fig. lip (a) is a finite size effect [22]). 

We corroborate these results with a calculation of the 
Chern number C . For a clean system: 

C = ^J^^ix{P{k){du,P{k\du,P{k)\}d^k, (1) 



where P(/c) is the /c-decomposition of the projector P 
onto the occupied states. In real-space: 

C = 27rz 5];^(0, Oi\P{- i\xx,P\ -z[f 2, P]] |0, a). (2) 

Eq. [2] is useful, as it allows one to treat finite disor- 
der. A classic result [23] states that the disorder average 
-27ri(^^(0,a|Pa;[[xi,Pa;],[f2,^a;]]|0,a))^, (Pa; = the 
projector onto the occupied states of H^) is an integer 
if Ep is in a region of localized EnS. This integer can 
change its value only if Ep crosses a region of extended 
states, a property that allows one to map the delocal- 
ized spectrum. Moreover, if R^ and Ep are chosen such 
that Ct^O at VK=0, then moving in any direction in the 
(Ep^ W^ plane from that initial point, C will eventually 
become zero for one of the reasons: 1) Pu;=0 if Ep is 
very negative, 2) Pu;=l if Ep is very positive or 3) all the 
states localize if W is too large. This implies the exis- 
tence of a region of extended states surrounding the CI 
phase and explains the phase diagram of Fig. [l] 

Unfortunately, Eq. [2] only makes sense in the thermo- 
dynamic limit. We derive a finite size real-space formula 
for C that converges exponentially fast to the thermody- 
namic limit. It does not involve twisted boundary condi- 
tions, which not only eliminates the problems associated 
with level crossings at Ep^ but allows us to compute C for 
large systems and many disordered configurations. Ad- 
ditionally, our formula for C requires only knowledge of 
the ground state. 

For clean CIs, C is computed using a discretized 
Brillouin-zone: /c-„,=^i^i+^2^2 5 ^1,2=1, • • • 7 ^7 l-^il = 
A = ^. The partial derivatives d]^^P are replaced by fi- 

nite differences SiP{kn)='^^CmP{kn -\- mAi) and the 
integration by a Riemann sum. Because the integrand 
in Eq. [l] is a periodic and analytic function, with an 
appropriate choice for the 5iP{kn) approximation, the 
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FIG. 2. (A) EtS and the variance of the level spacings (left panels), EnS (mid panel) and C number (right panels) for (A) 
C = 0.3i CI and (B) ^ = 0.3 normal insulator both with W—Z. The EtS and C were computed for seven ^f's as indicated. 
The blue/red data correspond to calculations on 30x30/40x40 lattices. C was identically zero for the normal insulator. The 
vertical dotted line in the EtS plots marks the 0.178 value. Panels (1-3) show the histograms of the level spacings collected at 
the entanglement energies marked by arrows. Overlayed, is the Wigner-Dyson distribution Pgue(s). 



discretized formula converges exponentially fast to the 
continuum limit. To obtain our real-space representa- 
tion, we note that the discretized C formula can be 
written as — zTr{P[(5iP, (52P]}, where the trace is over 
the whole Bloch basis \hnOL) and P is the full projector: 
P = ^ \hy^OL)P(y,^(hr\)l^y^^\. Slucc thc trace is invari- 
ant to a change of basis, we express this trace in the dual 
real-space basis |cc,a). The result is similar to that of 
Eq. [2]but with the substitution (c^ = — c_^): 

^ ^ — imA.iX -pAmA.iX 



Xij 1 
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CqJlC 
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If c^'s are chosen so that: x — ^J^^_jsj/2 
O(A^), the above substitution leads to exponentially 
small O(A^) errors. Together with the localization of P, 
this leads to an exponentially fast converging formula. 

For a clean CI (C=0.3i), the formula gives C= 
0.9999998/0.999999998 for a 30x30/40x40 lattice. The 
values at W=3 are shown in Fig. [2]A. These calculations 
were performed for the 30x30 and 40x40 lattices, 10^ 
configurations and for seven Ep values. The disorder av- 
eraged C values (for the 40x40 lattice) are 0.9999, 0.998, 
0.85, 0.53, 0.17, 0.01 and 0.0001 for Ef=^, 0.5, 1.0, 1.5, 
2.0, 2.5 and 3.0, respectively. These values indicate the 
existence of a delocalized spectral region between E=1.0 
and 2.0, in good agreement with Fig. lA. The topological 
character of the CI survives disorder. 

We now ask if the topological character of a CI has a 
clear signature in the EtS of the state. We compute the 
reduced density matrix by cutting a torus shaped sample 
in two equal parts A and B, and then tracing out B. 
The system is non-interacting, so we can use the single- 



particle EtS. This involves [24] finding the eigenvalues 
£^rn of the one-particle correlation function Cxy = {c\.Cy)^ 
where cc, y are lattice sites of the section that is not traced 
out^ and the expectation value is taken in the ground- 
state of the system. The reduced density matrix can be 
decomposed in normal modes with energies related to 
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= ilog(i 
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where ak 



are normal mode operators. The e^'s are "entanglement 
energies." In the clean limit, we perform translationally 
invariant cuts and plot ^^'s as function of momentum 
along the cut, as in FiglsFa). When Ep is in the bulk 
gap, the ^rn's are primarily concentrated around 0, 1 and 
have little dispersion. These are either bulk states deep 
in region A (^^'^1), or deep in region B (^m^O)- We call 
the difference between the levels at ^=1 and the ones at 
^=0 the entanglement bulk gap. For a trivial insulator, 
this is the whole story. For a nontrivial insulator as in Fig 
p[a), an entanglement mode localized on the cut crosses 
the entanglement bulk gap, much like an edge state in 
the EnS. As such, the EtS can differentiate between a 
topological and a trivial insulator, even though we are 
looking only at the bulk ground-state wavefunction for a 
system without boundaries. 

We now add disorder. Unfortunately, when Ep is at 
half filing, the number of the ^rn levels that are away from 
or 1 (i.e. the entanglement "edge spectrum") is very 
small, about 15 for a 30x30 lattice. As such, the large 
spacings between these levels render them in the clean 
regime of a random matrix. These "entanglement edge" 
levels exhibit strong level repulsion as can be clearly seen 
by the naked eye in FigpTb). The disorder mixing of 
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FIG. 3. a) Entanglement spectrum for a translationally in- 
variant CI plotted vs. momentum along the cut b) Level re- 
pulsion of the "edge" entanglement spectrum with disorder. 
Inset shows the full EtS. 



these levels is small (unless we go to high disorder where 
the CI is destroyed) and the level-spacing variance, al- 
though small, differs from 0.178. A computation with 
a larger 40x40 lattice shows a decrease of the variance. 
This is consistent with results in the almost clean limit of 
the EnS, where the disorder energy perturbation of each 
state is smaller than the mean level spacing [22]. These 
levels exhibit level repulsion and are delocalized. 

We now shift Ep and recompute the EtS and the vari- 
ance of its level spacings. The results are shown in Fig. |2] 
for W=3. As we move Ef towards the delocalized EnS, 
we notice that the level statistics of the EtS (in e^) ac- 
quires an increasingly flat region of level spacings dis- 
playing a variance of 0.178. As the Ef is moved up 
from half filling, the EtS becomes more and more dif- 
fusive (i.e. departs from the clean limit of the half-filled 
finite size problem). Delocalized bulk levels (which are at 
large negative or positive e) start moving in. When Ef 
sits right on top of the delocalized EnS, the whole EtS 
becomes delocalized and has variance extremely close to 
the Wigner-Dyson surmise of 0.178. The histograms of 
the level spacings collected from small windows at three 
widely spaced entanglement energies show well defined 
distributions matching closely the Pgve{s). The observa- 
tion of this delocalization plateau in the EtS correspond- 
ing to the ground-state of the system filled up to the ex- 
tended state energy is our main result. As Ef is moved 
above the extended states and into the region of the triv- 
ial Anderson insulator, the entanglement spectrum starts 
to become localized, with the spectrum near the center 
starting first. In contrast, for a trivial insulator, the EtS 
never has regions of level repulsion for any EtS of any 
ground-state. The calculations in Fig. |2] were performed 
for 30x30 and 40x40 lattices which give similar results 
although with reduced noise for the latter. 

The behavior of the EtS can be understood by noticing 
that Cxy is identical to the projection operator onto the 
occupied states, which can be regarded as a flat-band 
Hamiltonian. EtS is then identical to the spectrum of 
this flat-band Hamiltonian with open BC at the edge of 



the untraced region. The periodic BC fiat band Hamilto- 
nian has eigenvalues at and 1, and the open BC Hamil- 
tonian will have most of its eigenvalues close to or 1 
as well. If the original system was an insulator, then 
the eigenvalues between and 1 are separated by an en- 
tanglement gap, but if topologically nontrivial, the spec- 
trum also has edge modes crossing the entanglement gap. 
This is an inescapable conclusion that requires no calcu- 
lation. The projector onto occupied states satisfies the 
same symmetries as the original Hamiltonian. The pro- 
jector Hamiltonian has long-range hoppings when Ef is 
at the mobility edge, and we conjecture that this results 
in the flat, delocalized behavior of the full EtS. 

In conclusion, the EnS, EtS, and a new finite-size 
Chern number formula in the presence of strong disorder 
yield matching results and can be used to characterize 
the CI to Anderson-insulator transition. We found that 
all the levels of the EtS of a CI groundstate filled up to 
the edge of the mobility gap exhibit level repulsion con- 
sistent with the Wigner Dyson distribution (the many- 
body EtS matrices belong however to Wishart ensembles 
rather than Unitary ones). This delocalized plateau in 
the entanglement spectrum could be used to gain infor- 
mation about the many-body localization problem. 
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